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We search for hadronic decays of a light Higgs boson (A ) produced in radiative decays of an 
T(2S) or T(3S) meson, T —¥ "/A . The data have been recorded by the BABAR experiment at the 
T(3S) andT(25*) center of mass energies, and include (121. 3± 1.2) x 10 6 T(3S) and (98.3±0.9) x 10 6 
T(2S) mesons. No significant signal is observed. We set 90% confidence level upper limits on the 
product branching fractions B(T(nS) —¥ jA°) ■ B(A° — > hadrons) (n — 2 or 3) that range from 
1 x 10~ 6 for an ,4° mass of 0.3 GeV/c 2 to 8 x 10~ 5 at 7 GeV/c 2 . 

PACS numbers: 14.80.Da, 14.40.Pq, 13.20.Gd, 12.60.Fr, 12.15.Ji, 12.60.Jv 



A light CP-odd Higgs boson is expected in extensions 
to the Standard Model such as non-minimal Supersym- 
metry [l|. Light, in this context, means a mass less than 
that of the T(1S) meson. Such a Higgs boson could be 
produced in radiative decays of the T(nS) mesons 
Y(nS) —> 7A , where in this analysis, n = 2 or 3. BABAR 
has previously searched for this process where the A de- 
cays to muons [3] , taus [4] , or invisibly 0, || . CLEO has 
used its T(1S) data sample to search in the muon pair 
and tau pair final states [7[ . BABAR. has also searched for 
violations of lepton universality in T(1S) decay [8| , which 
could arise if the A has the expected quantum numbers 
J PC = CT+ and mixes with the r] b (lS) @. 

Supersymmetry models in which tan 2 /3 is not small 
predict that the A will decay predominantly into the 
heaviest kinematically available down-type fermion pair. 



The earlier experimental results have ruled out much of 
the parameter space 13, 11 1 . Regions not excluded tend 
to be dominated by hadronic decays, including decays to 
gluon pairs, gg, at smaller tan 2 (3, and to charm quark 
pairs, cc, at higher A mass. 

This analysis searches for hadronic decays of the A 
in the mass range 2m T < m A o < 7 GeV/c 2 without at- 
tempting to specify the underlying partons to which the 
A decays. The analysis nominally assumes that the A 
is CP-odd, but also relaxes this assumption to obtain 
results without specifying the CP state. 



The data were collected by the BABAR detector [12 1 
at the PEP-II asymmetric-energy e + e~ collider at the 
SLAC National Accelerator Laboratory. They consist 
of 27.9 fb _1 at the center of mass (cm.) energy of 
the Y(3S) and 13.6 fb -1 at the T(2S), correspond- 
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ing to N 3S = (121.3 ± 1.2) x 10 6 T(3S) and N 2S = 
(98.3 ± 0.9) x 10 6 T(25) mesons. We also use a con- 
tinuum (i.e., non-X(rtS')) background sample consisting 
of 78.3 fb _1 of data collected at the cm. energy of the 
T(45), plus 11.8 ftT 1 of data recorded 30-40 MeV be- 
low the r(2S), r(3S), or Y(4S) cm. energies. All of the 
data used here were recorded after the installation of an 
upgraded muon identification system [pj ]. 

Simulated signal events with various ^4° masses are 



used in the analysis. The EvtGen event generator [14 1 
is used to simulate particle decays. The A is simulated 
as a spin-0 particle, with equal branching fractions to 
whichever of gg, ss, and cc are kinematically available. 
Simulated events are produced both with and without 
the assumption that the A is CP-odd. JETSET [lj| is 
used to hadronize the partons, and Geant4 [iff is used 
to simulate the detector response. 

The search for the A uses hadronic events in which 
the full event energy is reconstructed. The selection cri- 
teria were optimized using simulated signal events and 
the continuum data set. The highest energy photon in 
each event is taken to be the radiative photon from the 
T(nS) decay. The ^4° candidate is constructed by adding 
the four-momenta of the remaining particles in the fol- 
lowing order. The first added are K® tt + tt~ candi- 
dates that have mass within 25 MeV/ c 2 of the true mass 
17], and whose reconstructed vertices are separated from 
the interaction point by at least three times the uncer- 
tainty on the vertex location. Charged hadron identifica- 
tion is then used to assign the proton, K , or mass 
to charged tracks. Tracks are labeled protons only if 
they are in the angular acceptance of the DIRC hadron 
identification system [la ], and if there is another track 
identified as an anti-proton. Neutral pion candidates 
are formed from pairs of photons, requiring the invari- 
ant mass of the photon pair to be between 100 MeV/c 2 
and 160 MeV/c 2 , and to have a 7r° energy greater than 
200 MeV. Finally, any remaining unused photons are 
added. Photons, including those used to reconstruct 
7T° mesons, are required to have a minimum energy of 
90 MeV. All energies and momenta are in the cm. frame. 

Events are required to have a radiative photon energy 
greater than 2.5 GeV (T(3S)) or 2.2 GcV (T(25)) and 
to have at least two charged tracks among the A decay 
products. The A mass resolution is improved by con- 
straining the radiative photon and all A decay products 
to come from a common vertex, and the sum of the pho- 
ton and j4° four-momenta to be that of the cm. system. 
To ensure that the full event energy is correctly recon- 
structed, the probability of the % 2 of the constrained fit 
is required to be greater than a value that ranges from 
at low A mass to 0.01 at m A a = 7 GeV/c 2 . Events 
in which m^o > 5 GeV/c 2 are rejected if the radiative 
photon, when combined with any other photon in the 
event, forms an invariant mass within 50 MeV/c 2 of the 
7T mass, or, for m A a > 6 GeV/c 2 , within 50 MeV/c 2 of 
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FIG. 1: Candidate mass spectrum in the (a) CP-all and (b) 
CP-odd analyses. The top curve in each plot is the on-peak 
data overlaid (in red) with the background fit described in the 
text, while the bottom curve (blue) is the scaled continuum 
data. The prominent initial state radiation resonances are 
labeled. 



the rj mass. 

Additional criteria are used to reject radiative Bhabha 
events, e + e~ — > je + e~ , or radiative muon pairs, e + e~ — > 
7/i + /i~. An event is rejected if it was identified as a 
Bhabha at the trigger level, if either of the two highest- 
momentum tracks is identified as an electron or a muon, 
or if the angle between the radiative photon and the 
second-highest momentum track is less than 1 radian. 
These criteria reject 96% of the continuum sample at a 
cost of 10-20% in signal efficiency, and, according to sim- 
ulation, reduce these backgrounds to negligible levels. 

The analysis proceeds along two parallel paths labeled 
"CP-all", in which no assumption is made on the CP 
nature of the A , and "CP-odd", in which it is assumed 
to be CP-odd. Events in which the A decays to ir + Tr~ 
or K + K~ are excluded from the CP-odd analysis. 

The analysis selects 371,740 events (CP-all) or 171,136 
events (CP-odd) in the combined T(2S) and Y(3S) ("on- 
peak") data set with 0.29 < m A o < 7.1 GeV/c 2 (Fig. [p. 

An A signal would appear as a narrow peak in the can- 
didate mass spectrum. The number of signal events at 
a particular hypothesis mass is computed as the number 
of events in a mass range ( "window" ) centered on that 
value, less the number of background events in the win- 
dow. The width of the window depends on the A mass 
resolution and was optimized along with the other selec- 
tion criteria. It varies for CP-all from 3 to 26 MeV/c 2 
as m A a increases from 0.29 to 7 GeV/c 2 . The CP-odd 
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windows are the same width as CP-all above 2 GeV/c 2 , 
but are larger at lower masses. 

Background events are from T(nS) decays and from 
continuum. Continuum, which is dominant, mostly con- 
sists of the initial state radiation (ISR) production of a 
light vector meson (clearly visible in Fig. [1} and non- 
resonant hadrons. The T(nS) backgrounds are primar- 
ily radiative decays to a light meson or non-resonant 
hadrons. At the highest ^4° candidate masses, there is 
an additional contribution from hadronic T(nS) decays 
in which a 7r° daughter is misidentified as the radiative 
photon. Simulation indicates that the fraction of BB 
events satisfying the selection criteria is negligible, so 
events recorded at the Y(AS) cm. energy can be used 
in the continuum sample. 

The number of background events is obtained from a 
fit to the data that contains three components: contin- 
uum, non-resonant T(nS) radiative decay, and resonant 
T(nS) radiative decay. The continuum component is the 
candidate mass spectrum of the continuum data set mul- 
tiplied by a normalization factor Cm (~ 0.5). Because 
the efficiency for detecting the ISR photon depends on 
cm. energy, Cm is not simply the ratio of integrated lu- 
minosities. It is left as a free parameter in the nominal 
fit, but, as described below, is fixed to a calculated value 
for systematic studies. The non-resonant T(nS) compo- 
nent is a 16-knot cubic spline, fixed to at the mini- 
mum .4° mass. The resonant component includes five 
rclativistic Breit-Wigner functions to represent the reso- 
nances for which CLEO saw some evidence in the study 
ofr(LS) -+jh+h- (h = novK) 18]: f (980), / 2 (1270), 
/£(1525), M1710), and / 4 (2050). The masses and widths 
are fixed 17J and possible interference between the res- 
onances is neglected in the fit. These resonances are all 
broad compared to an A signal. The spacing of the 
knots, typically 0.5 GeV/c 2 , is large enough that the cu- 
bic spline cannot conform to a narrow resonance. 

The background fit (Fig.[lJ has 21 free parameters, and 
is made to 1362 bins of width 5 MeV/c 2 , ranging from 
0.29 to 7.1 GeV/c 2 . The fit X 2 are 1268 (CP-all) and 
1293 (CP-odd) for 1341 degrees of freedom. Subtracting 
the normalized continuum mass spectrum from both the 
data and the fit gives the Y(nS) decay spectrum and the 
non-resonant and resonant radiative T decay components 
of the fit (Fig. HJ. 

The uncertainty on the background in each mass win- 
dow is both statistical and systematic. The systematic 
error is the sum in quadrature of the change in the total 
background arising from each of 17 alternative fits: the 
five nominal light resonances are removed one at a time, 
and eleven additional resonances are included one at a 
time. The eleven are established resonances 171 with 
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even total angular momentum, charge conjugation quan- 
tum number of +1, and Isospin 0. The seventeenth al- 
ternative fit is performed with C m fixed to the mid-point 
of the range of values found from four different methods 
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FIG. 2: A candidate mass spectrum after continuum sub- 
traction, overlaid with fit. (a) CP-all analysis; (b) CP-odd 
analysis. 



of determining it. Two of the methods are the nominal 
fits to the CP-odd and CP-all samples, and two use ISR- 
produced narrow resonances in four different final states: 
e + e~ — » yw, ijj — > it + TT~it a ; e + e~ — > 70, <fi ~^ K + K~\ 
e + e~ — > jJ/ip, J/tp — >> 4 charged tracks, with no ir°; 
and e + e~ — > 7 J/tp, J/ip ~ ^ 4 charged tracks, with one 
7T°. First, the number of each of these resonances is com- 
pared in on-peak and continuum data. Second, the same 
ratios are obtained using simulated samples of these ISR 
events, together with the calculated production cross sec- 
tions [19( and the recorded luminosities. The resulting 
value of Cm is 4.5% larger than nominal for CP-all, and 
2.7% for CP-odd. The fit qualities are good in all alter- 
native fits. The systematic errors are small compared to 
statistical errors except near resonances. 

The ^4° signal is evaluated at hypothesis masses that 
range from 0.291 GeV/c 2 to 7.000 GeV/c 2 in 1 MeV/c 2 
steps for the CP-all analysis (6710 mass hypotheses), and 
from 0.300 GeV/c 2 to 7.000 GeV/c 2 in 1 MeV/c 2 steps 
for CP-odd (6701 masses). Figure [3] shows the nomi- 
nal statistical significance of the resulting A signal, de- 
fined as the number of events divided by the statistical 
error, as a function of mass. The largest upwards fluc- 
tuations are 3.5<r at 3.107 GeV/c 2 for CP-all and 3.2cr 
at 0.772 GeV/c 2 for CP-odd. Including background sys- 
tematic errors, the significance of these two, which are 
located near the J/ip and p resonances respectively, are 
reduced to 2.8(7 and 2.2o\ The largest remaining fluctu- 
ations are 2.9cr at 1.295 GeV/c 2 for CP-all and 3.1er at 
4.727 GeV/c 2 for CP-odd. Figure 0] histograms the sta- 
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FIG. 3: Statistical significance (events divided by statistical 
error) of the ^4° signal as a function of mass, for (a) CP-all 
analysis, and (b) CP-odd analysis. 




A significance 

FIG. 4: Histogram of the statistical significance of the A 
signal for (a) the 6710 masses considered in the CP-all analy- 
sis, and for (b) the 6701 masses in the CP-odd analysis. The 
overlaid curve shows the distribution expected in the absence 
of signal. 



tistical significance of the signal measured at each mass, 
overlaid with the distribution expected in the absence of 
a signal. 

The signal extraction technique is studied using many 
simulated experiments. Each experiment consists of two 
candidate mass distributions, one for on-peak data, and 
the other for continuum. The continuum event distribu- 
tions are obtained from the full T(4S) data, which is 11 
times larger than the on-peak data set. The non-resonant 
T(nS) events are generated from a smooth threshold 



curve, and the resonant events are generated from rel- 
ativistic Breit-Wigner functions. The full signal extrac- 
tion is then performed. The average bias on the A signal 
yield is less than 1.5 events for all masses when there is 
no signal. 

These studies are also used to calculate the expected 
distribution of statistical significance in the absence of 
signal (Fig. 0} and to evaluate the significance of the 
largest apparent ^4° signals. The fraction of background- 
only CP-all simulated experiments that contain a fluctu- 
ation of nominal statistical significance > 3.5er is 33%. 
The fraction of CP-odd simulated experiments that con- 
tain a fluctuation > 3.2a is 63%. We therefore see no 
evidence of signal. The studies further indicate that 
large correlations between the resonant and non-resonant 
Y(nS) components make the uncertainties on the yields 
of the resonances unreliable. 

In the absence of a significant signal, we calculate a 
90% confidence level (CL) upper limit for each hypoth- 
esis mass on the product branching fractions B 3 s = 
B(Y(3S) -> jA°) ■ B{A° -> hadrons) and B 2S = 
B(T(2S) -> 7A ) • B(A° -> hadrons), assuming that the 
T(3S) and T(2S) decays are described by the same ma- 
trix element. This implies that B 2 s = B 3 s • ^3s/^ 2 s ■ 
R(rriA°), where T 3 $ and F 2 $ are the full widths of the 
T(3S) and T(2S) respectively, and R accounts for the 
difference in phase space. R is within a few percent of 
unity for all A masses. 

The calculation uses the relationship N = B + N^ s ■ 
B 3 s ' e> where N is the expected number of observed 
events for the given value of B%s, B is the expected 
background, N^ s = N 3S + N 2S ■ r 3S /r 2 s • R(m A o), 
and e is the signal efficiency. We calculate a likelihood 
C(Bss), defined as the probability of observing N or 
fewer events given that value of B 3 s, where N is the 
number actually observed. C(B 3 s) is obtained by inte- 
grating over the uncertainties in B, N 2 $, N3S, and e, 
which are assumed to be Gaussian. The 90% CL upper 
limit Bqq is calculated assuming a uniform prior above 0: 
j B9 " C(B 3S )dB 3S = 0.90 / °° C(B 3S )dB 3S . 

The efficiency is calculated using simulated events. 
The efficiency for the CP-all analysis ranges from a peak 
of 22% near m^o = 0.6 GeV/c 2 to less than 1% at high 
masses, while for the CP-odd analysis it ranges from 12% 
near 0.9 GeV/c 2 to less than 1% at high masses. 

The uncertainty on the efficiency is typically 11% (CP- 
all) or 7% (CP-odd) below the cc threshold, and 25% 
above. This includes contributions from uncertainty in 
tracking (1.5-3.5% depending on mass), photon and n° 
reconstruction (5-10%), and particle identification (3- 
5%), but the dominant contribution is due to the A 
decay branching fractions. This uncertainty is evaluated 
by varying the assumed branching fractions from 50% 
ss and 50% gg to 100% gg below the cc threshold, and 
from one-third each gg, ss, and cc to 50% cc and 25% 
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FIG. 5: 90% CL upper limits on product branching fractions 
(BF) (left axis) B(T(3S) -> jA°) ■ B(A° -> hadrons) and 
(right axis) B(T(2S) -> -yA°) ■ B(A° -> hadrons), for (a) CP- 
all analysis, and (b) CP-odd analysis. The overlaid curves in 
red are the limits expected from simulated experiments, while 
the blue curves are the limits from statistical errors only. The 
T(2S) limits do not include the phase space factor, which is 
at most a 3.5% correction. 

each of gg and ss above the cc threshold. The resulting 
systematic errors are 8% for CP-all or 4% for CP-odd 
below the cc threshold, and 21% above. The resulting 
90% CL upper limits are shown in Fig. [S] 

In conclusion, we have searched for hadronic final 
states of a light Higgs boson produced in radiative de- 
cays of the T(2S) or T(3S) and find no evidence of a 
signal. Upper limits on the product branching fraction 
B(T(nS) -> jA°) -B(A° -> hadrons) range from 1 x 1(T 6 
at 0.3 GeV/c 2 to 8 x 10~ 5 at 7 GeV/c 2 at the 90% CL. 
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